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We describe current- volt age characteristics I(V) of alkyl-ligated gold nanocrystals ^ 5 nm arrays 
in long screening length limit. Arrays with different alkyl ligand lengths have been prepared to tune 
the electronic tunnel coupling between the nanocrystals. For long ligands, electronic diffusion occurs 
through sequential tunneling and follows activated laws, as function of temperature a oc e~^°/^ and 
as function of electric field / oc e~^°/^. For better conducting arrays, i.e. with small ligands, the 
transport properties crossover to the cotunneling regime and follows Efros-Shklovskii laws as function 
of temperature a oc e~'^^^^^^^ ^ and as function of electric field / oc e~^^^^^^^ ^ . The data shows 
that electronic transport in nanocrystal arrays can be tuned from the sequential tunneling to the 
cotunneling regime by increasing the tunnel barrier transparency. 

PACS numbers: 73.23.Hk,73.21.La,71.30.+h 



Arrays of metallic, semiconducting, or magnetic 
nanocrystals (NC) with radii of 2 — 7 nm can now be 
synthesized [l|. Owing to their small self-capacitance, the 
charging energy for adding one electron per NC is large. 
Thus, these systems are ideally suited for the study of 
correlated electronic diffusion in presence of both disor- 
der and strong Coulomb interactions [I, [3|. 

In weakly conducting arrays, the large density of states 
at the Fermi level of metallic NC limits hopping to near- 
est neighbors and the large Coulomb energy opens a 
hard gap for electronic transport. These characteristic 
features led to an activated temperature dependence of 
the conductance where electron transport occurs by se- 
quential tunnel hopping [4]. While activated conductance 
was observed in granular metals [HI and in gold NC ar- 
rays at high temperature T > lOOi^^ [gI, H, the low tem- 
perature conductance of granular metals |8Ml2| and gold 
NC arrays [I, [l3l - [l7| has also been observed to follow the 
Efros-Shklovskii (ES) law : 



(1) 



This variable range hopping law is usually observed 
in doped semiconductors [18] where the wavefunction of 
dopant states decays exponentially and long range tun- 
neling between localized electrons is possible. 

The observation of the ES law in granular metals is 
puzzling as it implies long range tunneling of electrons 
beyond nearest neighbors. This issue was recently ad- 
dressed in a series of theoretical works [I, [l9l - [2H which 
established that effective long range hopping is possible 
because of cotunneling, where the electron hops to large 
distance through a string of virtual charge states. Such 
cotunneling phenomena is well known from the study of 
the transport properties of individual quantum dots [22]. 



In NC arrays, while the temperature dependence of the 
zero-bias conductance has been shown to follow Eq. [T] 
at low temperature, the question of the electric field de- 
pendence of electronic current in the cotunneling regime 
remains to be explored. 

In this letter, we present measurements of I(V) charac- 
teristics of gold NC arrays in the limit of long screening 
length, where the inter-electrode distance is small and the 
capacitive coupling of the NC with the gate is small. In 
this limit, the laws describing electron transport in the 
regime of cotunneling are known, they follow ES type 
formula both for the temperature dependence, Eq. [T]and 
the electric field dependence : 



(2) 



We study the evolution of the I(V) characteristics of 
the array as the electronic tunnel coupling between the 
NC is tuned by changing the ligand length. For small 
tunnel barrier transparency, the I(V) curves follow acti- 
vated laws as function of temperature and electric field; 
for large tunnel barrier transparency, the I(V) curves fol- 
low ES laws, Eq. [T]at low temperature and Eq. [2]at low 
electric field. 

Alkyl-ligated gold NC are synthesized by the digestive 
ripening method [j^]. This synthesis provides NC of ra- 
dius r :^ 2.5 nm with dispersion less than 5% in each 
batch. For this work, NC with alkylthiol-ligands of 6, 
8 and 12 carbon atoms have been prepared, so-named 
C^S^ CSS and C12S samples. The monodispersity of 
the NC allows for the preparation of hexagonal compact 
arrays with the Langmuir method [2^]. In a Langmuir 
trough filled with DI water, we first deposit a monolayer 
of dodecanethiol molecules. This amphiphile molecule 
improves the subsequent dispersion of gold NC. Then, a 
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FIG. 1. TEM image of dodecanethiol ligated gold NC self- 
organized in compact hexagonal array deposited on TEM 
grids by the Langmuir-Schaefer method. Panel a) Gold NC 
monolayer shown on a large scale (1 fim x 1 /jm). Inset) Gold 
electrodes on which the array is deposited. Panel b) Zoom on 
the hexagonal compact array. Panels c, d, e) one monolayer, 
two monolayers and three monolayers, respectively. 



compact array is obtained upon compression of the bar- 
riers. This process is monitored optically and through 
measurements of the surface pressure. 

After formation of the film, layers of NC are deposited 
with the Langmuir-Schaefer method; on substrates for 
conductivity measurements and on carbon grids for TEM 
characterization, shown Fig. [H Two layers are usually 
deposited on the substrate to reduce pin-holes in the ar- 
ray. We did not find any significant differences in the 
transport properties as the number of layers is changed 
from two to four layers. The substrates are p + + doped 
silicon covered with a silicon oxide layer 500 nm thick. 
On these substrates, we fabricated 2 /im width gold elec- 
trodes separated of the length L 150 nm. See in- 
set of Fig. [TJ CSS and C6S samples were measured as 
deposited on the substrate. C12S samples were dipped 
into isopropanol solution of alkyl-dithiol molecules during 
5 min to cross-link the NC and increase the tunnel trans- 
parency. I(V) measurements as function of temperature 
were performed after cross-linking in solutions of octane- 
dithiol, hexane-dithiol, butane-dithiol and ethane-dithiol 
molecules, respectively called C85'2, C65'2, C4S'2 and 
C2S2 samples. Despite the quality of the array, disorder 
remains. First, electrostatic disorder due to random off- 
set charges, second, tunnel barrier height disorder due to 
random ligands density and orientation. Both types of 
disorder are included within the theoretical model used 
below [21] to analyze the data, where electrostatic disor- 
der is modeled by grains with random potential and the 
tunnel matrix element is described by a random Gaussian 
variable. 

Figure [2] shows the I(V) curves of two as-deposited NC 
films CSS and C6S and one cross- linked film, C6S2. For 
all measured samples, the current increases non-linear ly 
with voltage, as typically observed for NC arrays in the 
regime of Coulomb blockade [3]. Figure [2li shows the tem- 
perature dependence of the zero-bias conductance of all 
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FIG. 2. Panels a) & b) I(V) curves of hexane-thiol (C6S) and 
octane-thiol (CSS) ligated NC arrays, respectively. Panel c) 
I(V) curve of dodecanethiol ligated NC crosslinked with hex- 
anedithiol ligands (C6S2). Panel d) Zero-bias conductance of 
the arrays with varying ligand lengths. The weakly conduct- 
ing arrays follow the activated law (dotted line). The more 
conducting arrays follow the ES law (dashed line). 



measured samples. Overall, the conductance of the ar- 
rays tends to increase as the ligands are shorter. This 
is expected as the conductance of N-alkanethiol ligands 
decreases exponentially fast with the number of carbon 
atoms TV as Gat - G'oe~^[25], where Go 2e^/h 
77.5 fiS. The zero-bias conductance of weakly conducting 
samples as CSS2 and CSS follows an activated temper- 
ature dependence aseq — (J^e~^^l^ from room tempera- 
ture to the lowest temperature (T 20 below which 
the conductance become smaller than the instrument sen- 
sitivity. The activation temperature Tq is similar for all 
weakly conducting samples, it is set by the Coulomb en- 
ergy Ec = e^/2Go where Go = ^.-Kse^r 0.8 aF is the 
self-capacitance of the NC where the dielectric coefficient 
e = 3 [2]. Using those parameters, we find that the activa- 
tion law with Tq = Ec/kB = 1110 can fit the zero-bias 
conductance of weakly conducting array as shown Fig. [2l 

In contrast, the conductance of best conducting ar- 
rays with shorter ligands, G4S'2, C2S2 and C6S^ deviates 
from this activation law at low temperature, T < SO K. 
This deviation cannot be due to a change in the acti- 
vation energy. Theoretical calculations [2l| give for the 
activation energy of arrays of dimensionless conductance 
g, Ea = Ec — {2gz/7r) * Ee/iln2, valid when gz « 1, 
where 2: = 6 is the coordination number and Eeh = 2Ec 
is the energy to create an electron-hole excitation in the 
system. As all the samples CnS2 are prepared from the 
same NC batch, there is no doubt that the NC size is 
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the same for all arrays, as are the self-capacitance and 
the Coulomb energy Eq- Furthermore, as the conduc- 
tivity of the arrays is very small, g = G/Gq ^ 10~^ at 
T :^ SOK for the best conducting arrays, this provides 
only a very weak correction to the activation energy and 
consequently Ea — Ec- 

We show now that this deviation can be understood 
as a consequence of inelastic cotunneling transport of 
electrons 1914211 26| . In this regime, the conductance 
should follow Eq. [H We find that this last formula fits the 
data with Tes = e^/{ATTeeQ < ^ > ks) ^ 8000 K. The 
average value of the localization length < ^ >= 0.1 6r 
is obtained from ^ = [E'l^/w^ikBTyg] ' ^^^^^ ^^^^ ^ 



1 [2l| . This last formula shows that the localization length 
changes slowly with temperature and conductance : from 
^ 0.23r {T = SO K, g = 2 X lO'^ ) to ^ 0.09r 
{T = 5 K, g = 4x lO'^ ). 

We calculate that the hopping distance rhop = 
V'e2^/(8A:BT7r££o), which increases with decreasing tem- 
perature, exceeds two NC below T ^ 50 K. This analysis 
confirms that electron transport occurs through cotun- 
neling below that temperature and is in agreement with 
the conclusions of previous works by Tran et al. on sim- 
ilar gold NC arrays [2j. We show now that an analysis of 
the electric field dependence of I(V) curves at low tem- 
perature allows to establish the signature of cotunneling 
transport in the array. 
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FIG. 3. Panels a) to c) show that the I(V) curves are self- 
similar, i.e. the curves appear the same on different scales. 
This implies that the apparent voltage threshold Vt indi- 
cated by arrows changes with the scale. This shows that 
Vt is not a relevant parameter for the transport properties 
of the array. In contrast, the Arrhenius plot - panel d) - 
shows that the electronic current follows an activated behav- 
ior I{V) oc e"^o/^ (dotted line). 



In previous works on gold NC arravs|2 Ul4l - ll6| . the non- 
equilibrium transport properties have been analyzed in 
the context of Middleton and Wingreen (MW) model (27| 
of percolating currents. One essential parameter of this 
model is the screening length A/r = Cq/Cg, where Co = 



Aneeor'^/d is the capacitive coupling of the NC with the 
gate electrode. Short screening lengths, A/r oc d/r 
0, are obtained in arrays strongly coupled to the gate 
electrode such as lithographically defined quantum dots 
[iij. In this limit, the MW model predicts the existence 
of a voltage threesold Vt for electron transport and a 
power law dependence of the current with voltage, / cx 
(V — Vt)^. As the effects of cotunneling on this model 
are not known theoretically Q, this prevents a reliable 
analysis of the effect of cotunnehng on I(V) curves. 

In our sample configuration, the screening length A ^ 
d 500 nm is larger than the electrode separation 
L = 150 nm, implying that the MW model does not 
apply fisl. In this long screening length regime, the laws 
describing the electric field dependence of the electronic 
current in the sequential and cotunneling regimes have 
been calculated theoretically [Hi |2l|. They mirror the 
temperature dependence and can be obtained by replac- 
ing the thermal energy ~ ksThy the electrostatic energy 

Experimentally, while an activated electric field depen- 
dence has been found for the conductance of granular 
metals [30], there is no report of the observation of the 
ES law Eq. [2]for the electric field dependence. 

Figure [3] shows data for a OSS' sample. On left pan- 
els, the data is shown on a linear plot at three different 
scales of increasing magnification where arrows indicate 
for every I(V) curve the apparent voltage threshold Vt. 
It clearly appears that Vt decreases upon increasing the 
plot magnification, which shows that Vt is not a relevant 
physical parameter. On the right panel, the same data 
on an Arrhenius plot shows that the current follows the 
activated law for sequential tunneling, Iseq — l^e~^^l^ 
where = kBT^/el is set by the Coulomb energy, / is 
the inter-nanocrystal distance on which the electron is 
submitted to the electric field, £ = cV/L^ where c takes 
into account the screening effects that reduce the actual 
electrical field from the nominal value V/L. A good fit of 
the data can be obtained with c = 0.15 and / = 1 nm, giv- 
ing for the activation electric field So = 1.3 x 10^ Vm~^. 
While there is a large uncertainty on the actual electric 
field applied on the electron because of the uncertainty 
on inter-nanocrystal distance / and on screening effects, 
we will see now that using the same values for the pa- 
rameters So and c, the ES law can fit the data in the 
cotunneling regime. 

Figure H] shows the I(V) curves for all samples. The 
I(V) curves of weakly conducting samples, C8S2 and 
C8S, follow the activated electric field dependence just 
described on the whole measured voltage range. In con- 
trast, the I(V) curves of best conducting arrays, C45'2, 
C2S2 and CG^, deviate from this law at low electric field. 
It follows the ES law Eql2]with the same characteristic 
electric field Ses = fo = 1-3 x 10^ Vm~^. Upon increas- 
ing the voltage, a crossover is observed at an electric field 
S '^10^ Vm~^ and the I(V) curve recovers the activation 



4 



Bias (V) 




1E-5 



1.5x10"^ 3.0x10"^ 4.5x10"^ 

1/£(\A^ m) 

FIG. 4. Panels a) shows the evolution of the I(V) curves of 
the arrays as the tunnel transparency between NC increases. 
For weakly conducting NC arrays, i.e. small tunnel barrier 
transparency, the I(V) curve follows an activated law, dotted 
line. As the conductance of the arrays increase, i.e. for larger 
tunnel barrier transparency, the I(V) curve follows an ES-type 
law, dashed line. The arrow indicates the field scale where 
the crossover between the two regimes occurs. Panel b) The 
thin red line is the hopping distance. The dotted line is the 
activated law, the dashed line is the ES law. As the hopping 
distance increases at low electric field, a crossover occurs from 
the activated law to ES law, as shown schematically by the 
thin blue line. 



law at high electric field. This is a remarkable observa- 
tion. We find that the electric field dependence mirrors 
the temperature dependence and that a single character- 
istic electric field value, set by the Coulomb energy, can 
describe the I(V) curves for a set of NC arrays whose 
room temperature conductance changes by three orders 
of magnitude. 

The same figure shows that the calculated hopping dis- 
tance, rhop = ^/e/{S^7^eeo)[21] increases rapidly about 
the electric field value where the crossover between the 
activated law and the ES law is observed. This con- 
firms our analysis, the crossover occurs between the co- 
tunneling regime at low electric field, where the hopping 
distance rhop ^ 5 x r at zero-bias, and the sequential 
tunneling regime, where rhop — r at high bias. 

To summarize, we studied the I(V) curves of gold NC 
arrays with long screening length. In this regime, the 
eflFects of cotunneling on the I(V) curves are known the- 



oretically, whose predictions can be compared to the 
data. As a function of temperature, the conductance 
follows an activated law at high temperature in the se- 
quential tunneling regime and an ES law at low tem- 
perature in the cotunneling regime. We find that the 
electric field dependence mirrors this temperature depen- 
dence. The I(V) curves follow an activated field depen- 
dence at high electric field in the sequential tunneling 
regime and an ES law at low electric field in the cotun- 
neling regime. The crossover from the activated law to 
the ES law as array conductance is increased can be un- 
derstood from the early calculations for linear NC chains 
where it was shown 31. [32j that the cotunneling current 
across j junctions increases with the junction conduc- 
tance G as / (X [■§^y ' 
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